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ABSTRACT 


The  lift,  drag,  and  inertia  coefficients  have  been 
determined  experimentally  for  various  tube  bundles  and  for 
two  cylinders  at  various  relative  positions  subjected  to 
harmonically  oscillating  flow.  The  force  coefficients  for 
the  inline  force  have  been  analyzed  through  the  use  of 
Morison's  Equation  and  Fourier  Analysis.  The  transverse 
force  or  the  lift  force  has  been  expressed  in  terms  of  its 
maximum  value . 

The  results  have  shown  that  the  interference  between  the 
cylinders  can  give  rise  to  complex  flow  patterns  and  to  un¬ 
expectedly  large  force  transfer  coefficients.  The  results 
with  two-cylinder  experiments  have  shed  considerable  light 
on  the  wake  interference  in  harmonic  flows. 
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I .  INTRODUCTION 


A  body's  resistance  to  flow  is  strongly  affected  by  what 
surrounds  it.  When  two  bodies  are  in  close  proximity,  not 
only  the  flow  about  the  downstream  body  but  also  that  about 
the  upstream  body  may  be  influenced.  Examples  include  con¬ 
denser  tubes  in  heat  transfer;  columns  in  pressure  suppression 
pools  of  nuclear  reactors;  risers,  piles,  and  other  tubular 
structures  in  off-shore  engineering;  turbine  and  compressor 
blades  in  mechanical  or  aerospace  engineering;  and  high-rise 
buildings,  cooling  towers,  and  transmission  lines  in  civil 
engineering.  The  quantification  of  the  interference  effects 
in  terms  of  the  pressure  distribution,  lift  and  drag  forces 
on  individual  members,  vortex  shedding  frequency,  and  the 
dynamic  response  of  members  of  the  array  in  terms  of  the 
governing  flow  and  structural  parameters  constitute  the 
essence  of  the  problem. 

There  are  infinite  numbers  of  possible  arrangements  of 
two  or  more  bodies  positioned  at  right  or  oblique  angles  to 
the  approaching  flow  direction.  In  wavy  or  time-dependent 
flows  one  needs  the  lift,  drag  and  inertia  coefficients  for 
all  members  of  the  array.  Evidently,  the  members  of  the 
array  may  not  be  all  parallel  and  normal  to  the  flow.  The 
quantification  of  the  flow  interference  on  lift,  drag  and 
inertia  coefficients  for  cylinders  with  relative  inclinations 
and  spacings  in  a  design-wave  environment  is  an  exceedingly 
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complex  problem.  In  the  absence  of  data  on  the  interaction 
between  drag  and  inertia  coefficients  in  wavy  and  harmonic 
flows  for  cylinder  groups  one  is  tempted  to  use  steady-flow 
results  for  the  drag  coefficient  and  the  unseparated  potential 
flow  results  for  the  inertia  coefficient.  The  inertia  coef¬ 
ficient  for  a  group  of  cylinders  in  inviscid  unseparated 
flow  may  be  obtained  through  the  use  of  the  method  of  images 
/  1,  2_/  or  through  the  use  of  the  linear  potential  theory 
including  wave  diffraction  /~ 3 ,  4_/,  These  analyses  do  not 
deal  with  the  effects  of  separation  and  vortex  shedding.  Con¬ 
sequently,  the  results  are  more  appropriate  to  the  determina¬ 
tion  of  earthquake  forces  and  wave  forces  on  large  bodies 
rather  than  to  the  evaluation  of  the  inertial  components  of 
the  force  in  the  drag/inertia  dominated  regime. 

It  is  evident  from  the  foregoing  that  experiments  must  be 
carried  out  with  carefully  selected  tube  bundles  and  cylin¬ 
der  combinations  in  order  to  develop  some  understanding  of 
the  flow  interference  and  to  provide  data  for  body  combina¬ 
tions  of  practical  importance. 
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II.  REVIEW  OF  PREVIOUS  INVESTIGATIONS 


A  careful  review  of  flow  interference  between  two  circular 
cylinders  in  various  arrangements  in  steady  flow  has  been  pre¬ 
sented  by  Zdravkovich  /  5  /  where  an  extensive  list  of  re¬ 
ferences  may  be  found.  Numerous  studies  have  shown  that  the 
changes  in  drag,  lift  and  vortex  shedding  are  not  necessarily 
continuous.  In  fact,  the  occurrence  of  a  fairly  abrupt  change 
in  one  or  all  flow  characteristics  at  a  critical  spacing  is 
one  of  the  fundamental  observations  of  flow  interference  in 
cylinder  arrays. 

For  the  tandem  arrangement  (one  cylinder  behind  the  other, 
see  Figure  1) ,  it  has  been  shown  that  at  relative  spacings 
S/D<3.5  there  is  a  strong  mutual  interference  between  the  two 
cylinders.  This  critical  distance  increases  with  the  bluff ness 
of  the  bodies.  For  two  plates  in  tandem,  the  critical  spacings 
is  about  10  times  the  plate  width  /  6_/.  In  general,  the 
tandem  arrangement  has  a  strong  effect  on  the  drag  and  is  sen¬ 
sitive  to  spacing.  The  upstream  cylinder  takes  the  brunt  of 


Fig.  1  Tandem  arrangement  of  two  cylinders 


the  burden  and  the  total  drag  for  the  group  is  smaller  than 
the  sum  of  the  drag  forces  acting  on  each  cylinder  in  isola¬ 
tion  in  a  tunnel  with  the  same  blockage.  The  drag  coefficient 
for  the  downstream  cylinders  is  reduced  partly  by  shielding 
and  partly  by  the  occurrence  of  earlier  transition  in  the 
boundary  layers  due  to  the  turi ulence  in  the  approaching  flow. 


D  h*- 


Fig.  2  Side-by-side  arrangement  of  two  cylinders 

The  side-by-side  arrangement  (see  Figure  2)  exhibits  a 
discontinuity  in  the  flow  and  measured  forces  for  spacing 
ratios  smaller  than  a  critical  value  (S/D<2.2  for  two  circu¬ 
lar  cylinders  and  about  4  two  plates) .  The  observed  dis¬ 
continuity  (a  switching  phenomenon  and  bistable  nature  of 
the  flow  between  plates)  results  in  two  values  of  the  drag 
coefficient.  This  results  from  the  mutual  interference  of 
the  vortices  on  the  adjacent  sides  of  the  vortex  streets. 

The  bistable  nature  of  the  flow  may  be  avoided  by  increasing 
the  size  of  one  of  the  cylinders.  The  sum  of  the  bistable 
high  and  low  drag  is  often  less  than  twice  the  drag  of  the 
single  cylinder.  For  two  plates  in  side-by-side  arrangement, 
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the  said  sum  may  be  about  10  percent  larger  than  twice  the 
drag  of  the  single  plate.  This  increase  strongly  depends  on 
the  width  of  the  flow  field  relative  to  the  total  width  of  < 

the  bodies,  i.e.,  on  the  blockage  ratio.  For  relatively 
large  blockage,  the  flow  is  forced  through  the  openings 
between  the  cylinders  rather  than  around  the  total  configura¬ 
tion.  Consequently,  the  use  of  ordinary  shielding  and  block¬ 
age  factors  for  groups  of  cylinders  is  meaningless.  Experi¬ 
ments  must  be  conducted  in  channels  or  tunnels  with  very  small 
blockage  in  order  to  obtain  valid  force-transfer  coefficients. 

There  is  at  present  no  means  to  separate  the  flow  interference 
effects  from  the  blockage  effects  for  groups  of  cylinders. 

Additional  work  with  tube  arrays  in  steady  flows  has  been 
reported  by  a  number  of  people  /~7  -  14 Ross  /~ 15  7  con¬ 
ducted  large  scale  wave  tank  tests  for  the  case  of  one  cylin¬ 
der  on  each  side  of  the  test  cylinder  in  the  range  of  critical 
Reynolds  numbers.  His  results  appear  to  indicate  that  the 
wave  force  increases  significantly  only  when  the  spacing 
between  two  cylinders  is  less  than  about  one  diameter. 

The  foregoing  studies,  conducted  mostly  in  steady  flows, 
show  that  the  sum  of  the  individual  drag  forces  for  isolated 
cylinders  is  often  greater  than  the  total  force  acting  on  the 
array  as  long  as  the  local  boundaries  do  not  constrain  the 
flow.  They  further  show  that  the  results  obtained  from  the 
tandem  and  side-by-side  arrangements  cannot  be  generalized  to 
predict  the  combined  interference  resulting  from  the  staggered 
arrangement. 
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Gibson  and  Wang  /  16_/  carried  out  two  different  experi¬ 
ments  to  determine  the  added  mass  of  a  series  of  tube  bundles. 
The  bundles  consisted  of  tubes  of  uniform  diameter  d  arranged 
either  in  a  square  configuration  or  a  circular  configuration. 

In  the  first  series  of  experiments,  they  towed  the  model  of  a 
pile  cluster  under  linear  acceleration.  In  the  second  series, 
they  have  vibrated  the  model  at  its  own  natural  frequency.  For 
both  cases,  they  have  calculated  the  added  mass  through  the 
use  of  the  measured  force  and  acceleration  and  plotted  them  as 
a  function  of  the  'solidification  ratio'  defined  by  Zd^/irD 
where  D  is  the  pitch  diameter  of  the  bundle.  Their  results 
have  shown  that  the  added  mass  increases  sharply  after  the 
solidification  ratio  reaches  the  value  of  0.4  to  0.5.  Beyond 
this  value,  the  volume  enclosed  rather  than  the  volume  dis¬ 
placed  by  the  structure  becomes  important.  This  result  is 
disputed  /  17_/  and  the  results  of  both  series  of  tests  are 
no  more  applicable  to  separated  wavy  or  oscillatary  flows 
about  tube  bundles  than  those  predicted  from  the  potential 
theory  with  or  without  diffraction  effects. 

Relatively  few  studies  have  been  carried  out  with  oscil¬ 
lating  tube  bundles  /— 18,  19_7.  Bushnell  /  19_/  determined 
the  interference  effects  on  the  drag  and  transverse  forces 
acting  on  a  single  member  of  a  two-cylinder  configuration 
through  the  use  of  a  pulsating  water  tunnel.  He  did  not 
evaluate  the  drag  and  inertia  coefficients  through  the  use 
of  a  suitable  method,  e.g.  Fourier  averaging.  Instead,  he 
picked  out  the  maximum  force  values  which  occurred  in  each 
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half  cycle  and  averaged  them  over  ten  consecutive  values  so  as 
to  obtain  a  mean  maximum  force  for  each  flow  direction.  The 
results  have  shown  that  the  presence  of  neighboring  cylinders 
significantly  affects  the  forces  on  an  individual  cylinder  of 
an  array  and  the  interference  effect  increases  with  increasing 
relative  flow  displacement.  The  maximum  drag  force  on  shielded 
cylinders  was  reduced  relative  to  an  exposed  cylinder  by  up  to 
50  percent.  Bushnell  has  suggested  that  a  design  using  a  high 
Reynolds  number  single-cylinder  drag  coefficient  applied  through¬ 
out  the  array  would  have  an  extra  margin  of  safety  against  maxi¬ 
mum  drag  loading  due  to  interference  effects.  The  transverse 
force  could  be  3  to  4  times  larger  for  interior  array  positions 
than  that  for  a  single  cylinder.  Thus,  a  cylinder  array,  such 
as  a  riser,  supported  at  regular  intervals  may  exhibit  very 
complex  dynamic  behavior .  Some  of  the  members  between  supports 
may  undergo  violent  transverse  oscillations. 

It  is  evident  from  the  foregoing  that  there  is  very  little 
systematic  data  for  tube  bundles  and  two  cylinder  interference. 

In  view  of  this  fact  is  was  decided  to  carry  out  experiments 
with  two  bundle  configuration  consisting  of  12  and  9  outer 
tubes,  respectively.  Subsequent  to  this  investigation  a 
detailed  study  of  the  interference  effects  of  two  cylinders 
in  harmonic  flow  was  undertaken. 


III.  EXPERIMENTAL  EQUIPMENT  AND  PROCEDURES 


A.  OSCILLATING  FLOW  TUNNEL 

Experiments  were  carried  out  in  a  large  U-shaped  oscillating 
water  tunnel  (see  Figure  3) .  It  has  been  used  extensively  at 
this  institution  over  the  past  6  years  /— 20 _J .  Only  salient 
features,  most  recent  modifications,  as  well  as  the  adaptation 
for  this  work,  are  briefly  described  here. 

The  length  of  the  tunnel  has  been  increased  from  30  feet  to 
35  feet  and  its  height  from  16  feet  to  24  feet.  Previously,  a 
butterfly-valve  arrangement  at  the  top  of  one  of  the  legs  was 
used  to  initiate  the  oscillations.  In  1978,  the  tunnel  was 
modified  so  that  the  oscillations  could  be  generated  and 
maintained  indefinitely  at  the  desired  amplitude.  For  this 
purpose  the  output  of  a  2  HP  fan  was  connected  to  the  top  of 
one  of  the  legs  of  the  tunnel  with  a  large  pipe.  A  small 
butterfly  valve,  placed  in  a  special  housing  between  the  top 
of  the  tunnel  and  the  supply  line,  oscillated  harmonically 
at  a  frequency  equal  to  the  natural  frequency  of  the  oscilla¬ 
tions  in  the  tunnel.  The  oscillation  of  the  valve  was  per¬ 
fectly  synchronized  with  that  of  the  flow  through  the  use  of 
a  feedback  control  system.  The  output  of  a  pressure  trans¬ 
ducer  (sensing  the  instantaneous  acceleration  of  the  flow) 
was  connected  to  an  electronic  speed-control  unit  coupled  to 
a  DC  motor  oscillating  the  valve  plate.  The  circuit  main¬ 
tained  the  period  of  oscillations  of  the  valve  within  0.001 
seconds.  The  fluid  oscillated  with  a  period  of  6.000  seconds. 


The  amplitude  of  oscillations  was  varied  by  constricting  or 
enlarging  an  orifice  at  the  exit  of  the  fan.  The  flow  oscil¬ 
lated  at  a  given  amplitude  as  long  as  desired. 

Two  identical  force  transducers,  one  at  each  end  of  the 
cylinder,  were  used  to  measure  the  instantaneous  in-line  and 
transverse  forces.  A  special  housing  was  built  for  each  gage 
so  that  it  can  be  mounted  on  the  tunnel  window  and  rotated  to 
measure  either  the  in-line  or  the  transverse  force  alone. 

The  calibration  of  the  each  gage  was  accomplished  by 
hanging  loads  in  the  middle  of  the  cylinder  after  setting  both 
gages  to  sense  only  the  transverse  force  (here  in  the  vertical 
direction) .  Repeated  calibrations  have  shown  that  (a)  the 
gages  were  absolutely  linear  up  to  500  pounds;  (b)  the  gages 
yielded  the  same  signal  for  loads  applied  either  upward  or 
downward;  and  that  (c)  the  gages,  together  with  the  electronic 
system  to  which  they  are  attached,  were  capable  of  sensing 
loads  as  small  as  0.02  pounds. 

The  displacement,  velocity,  and  acceleration  of  the  fluid 
in  the  tunnel  were  obtained  through  the  use  of  several  methods 
described  in  /  20 7* 

The  data  were  recorded  in  analog  form  and  then  digitized 
and  processed  through  the  use  of  an  HP9845  computer-digitizer 
system. 

B.  TEST  BODIES 

The  tube  bundles  were  made  of  aluminum  rods  and  mounted 
between  two  circular  stainless  steel  plates.  Two  configura¬ 
tions  were  chosen  for  test  purposes.  These  configurations 
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Fig.  4b  Riser  2  configuration 
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roughly  correspond  to  the  most  commonly  used  riser  shapes  in 
the  off-shore  industry.  Figures  4a  and  4b  show  the  two  riser 
configurations  together  with  their  characteristic  dimensions. 

The  two  cylinder  configuration  was  arranged  in  a  similar 
manner  (see  Figure  5) .  The  active  cylinder  was  mounted  on  the 


Fig.  5  General  arrangement  of  two  cylinders 


force  transducer  as  described  previously.  The  passive  cylin¬ 
der  spanned  the  entire  test  section  and  was  rigidly  mounted 
on  the  tunnel  walls.  Both  cylinders  were  smooth  and  3  inches 
in  diameter.  The  angle  a  was  varied  from  zero  to  90  degrees 
in  30  degree  increments.  The  N  in  cylinder  spacing  ND  was 
varied  from  1.5  to  3.5.  The  in-line  and  transverse  forces 
acting  on  the  active  cylinder  were  measured  and  analyzed  in 
terms  of  the  appropriate  lift  coefficients  and  the  drag  and 
inertia  coefficients  through  the  use  of  the  Fourier  analysis. 
The  details  of  the  definition  of  the  various  force  transfer 
coefficients  are  described  in  the  following. 
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IV.  FORCE  TRANSFER  COEFFICIENTS 


A.  TUBE  BUNDLES 

The  in-line  force,  based  on  Morison's  equation  is  written 
as 


2 

F  =  -0.5/C.U  L  Z  D.  !  Coscjt  j  Coscat 

dm  1  1  1 

(1) 

+  0.25tt/LC  Z  D2U  jSin^t 
m  i  m 


in  which  U  is  the  maximum  velocity  in  a  given  cycle. 
Denoting 


D  =  Z  D. 
a  i 


and  inserting  in  equation  (1) ,  one  has 
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in  which  T  =  2tt/j  . 

For  an  oscillating  flow  represented  by  U  =  -U  cos^t  the 
averages  of  C^  and  Cm  are  given  by  /~ 20_/ 
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in  which  F  represents  the  measured  force, 
m 

The  lift  coefficient  is  defined  as 


C 


L 


0.5/LD  U2 
e  m 


(5) 


in  which  F^  represents  the  amplitude  of  the  transverse  force 
and  L,  the  length  of  the  cylinder. 

The  vortex  shedding  frequency  for  the  entire  bundle  was 
expressed  in  terms  of  the  Strouhal  number  defined  by 


St 


f  D 
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in  which  f^  represents  the  frequency  of  vortex  shedding;  T, 
the  period  of  flow  oscillation;  and  K,  the  Keulegan-Carpenter 
number  defined  by 


U  T 
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B.  TWO  CYLINDER  ARRANGEMENT 

The  in-line  force  acting  on  the  active  cylinder  is  not  in¬ 
dependent  of  the  direction  of  the  fluid  motion  because  of  the 
interference  effect  of  the  passive  cylinder.  In  the  interval 
tt/2 <ujt<3iT/2  the  direction  of  flow  is  opposite  to  that  in  the 
interval  3tt/2  <cut<5iT/2.  With  this  fact  in  mind  Morison's  equa¬ 
tion  was  written  as 

2 

F  it  D  1 
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0.5/D  U  m  d 

m 


and  C  and  C,  were  determined  from 
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in  which  and  represent  the  inertia  and  drag  coefficients 

for  the  interval  tt/2  <oit <3tt/2  ;  and  C  ^  and  /  the  corresponding 
coefficients  for  the  interval  3tt/2  <ajt  <5tt/2  .  In  general  these 
two  sets  of  in-line  force  coefficients  are  not  identical  accept 
when. a=90°. 

The  evaluation  of  the  transverse-force  coefficients  re¬ 
quires  a  distinction  not  only  in  the  direction  of  flow  but 
also  in  the  direction  of  the  transverse  force.  As  the  flow 
proceeds  in  a  given  direction  the  lift  force  acting  on  the 
active  cylinder  does  not  remain  symmetrical  primarily  because 
of  the  proximity  of  the  passive  cylinder,  situated  asymmetri¬ 
cally  relative  to  the  direction  of  the  flow.  It  is  also  re¬ 
cognized  that  the  random  nature  of  the  vortex  shedding  does 
not  allow  a  perfectly  symmetric  lift  force  even  in  the  absence 
of  any  interference  effects.  Consequently,  one  must  define 
four  lift  coefficients.  In  the  interval  tt  /  2  <  cot.  <  3  tt/ 2  ,  one 
has  (see  Fig.  6a) 
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L1  0.5/U2  DL 
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and  in  the  interval  3it/2 <oit <5^/2 ,  (see  Fig.  6b) 
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and 


where 
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Fl1  through  F  are  shown  in  Figures  6a  and  6b. 


Passive  Cyl.-^ 


Flo1 


o.  > 
*  — 


4  °Ll 


Active  Cyl.  - 1 


r  i 


'L2 


Fig.  6a  Definition  sketch  for  lift  forces 
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Fig.  6b  Definition  sketch  for  lift  forces 
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The  Strouhal  number  was  defined  in  a  manner  similar  to  that 


previously  described  and  evaluated  through  the  use  of 


in  which  i  =  1,  ..  .4,  indicating  the  lift  force  directions 

shown  in  Figures  6a  and  6b. 

C.  GOVERNING  PARAMETERS 

It  has  been  shown  previously  /  20 _7  that  the  Fourier- 
averaged  lift,  drag,  inertia  coefficients  for  a  smooth  cylin¬ 
der  immersed  in  a  harmonically  oscillating  flow  depend  on 
the  Reynolds  number  and  the  Keulegan-Carpenter  number.  In 
other  words 


(C , ,  C  ,  C  } 
d  m  L 


U  D  U  T 
m  m 

v  '  D 


In  terms  of  a  frequency  parameter  3  defined  by 
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one  has 


tCd'Cm'CL>  -  Fi  <§T  ' 


For  the  two  cylinder  arrangement  one  has  to  introduce  two 
additional  parameters  namely,  the  angle  a  of  the  line  joining 
the  center  of  the  two  cylinders  and  the  relative  distance 
between  the  cylinders.  Thus,  all  the  force-transfer  coeffi¬ 
cients  cited  previously  depend  on  the  Keulegan-Carpenter 
number,  Reynolds  number  (or  the  frequency  parameter) ,  the 
angle  a,  and  the  relative  distance  factor  N  (see  Figure  5) . 
Thus  one  has 


C^(a  coefficient)  =  F^(K,B,a,N) 


28 


V.  DISCUSSION  OF  RESULTS 


A.  TUBE  BUNDLES 

The  drag  and  inertia  coefficients  for  the  two  riser  con¬ 
figurations  are  shown  in  Figures  7a  and  7b  as  a  function  of 
the  Keulegan-Carpenter  number  for  constant  values  of  g.  The 
drag  coefficient  decreases  gradually  with  increasing  K.  The 
inertia  coefficient  increases  with  increasing  K  and  reaches 
values  as  large  as  6.  The  use  of  the  potential  flow  theory 
for  unseparated  flow  yielded  /~21_/  an  average  inertia 
coefficient  C*  =2.17  for  the  two  configurations.  C*  was 
defined  by 


C*  =  (IC  ,D2)/ZD2 
m  mi  i  i 


The  comparison  of  C*m  with  that  obtained  experimentally  has 
shown  that  some  fluid  mass  is  entrapped  within  the  array  and 
that  neither  the  potential  theory  nor  the  diffraction  analy¬ 
sis  can  adequately  describe  the  behavior  of  the  complex 
separated  flow  through  the  bundle. 

A  comparison  of  the  Figures  7a  and  7b  shows  that  the  drag 
and  inertia  coefficents  as  defined  here  are  nearly  indepen¬ 
dent  of  the  particular  arrangement  and  the  number  of  tubes 
comprising  the  bundle. 
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Experiments  with  different  tube  arrangements  /  21_/  have 
shown  that  the  force  coefficients  are  independent  of  the 
Reynolds  number  within  the  range  of  Reynolds  and  Keulegan- 
Carpenter  numbers  reported  herein. 

The  reason  for  the  dependence  of  and  on  K  is  thought 
to  be  the  dependence  of  the  interaction  of  the  wakes  of  the 
outer  and  inner  pipes.  The  vortices  in  the  wake  of  a  given 
cylinder  loose  about  70  percent  of  their  strength  within  10 
cylinder  diameters.  Thus,  for  small  values  of  K,  the  vortices 
generated  by  a  small  tube  at  the  center  front  of  the  bundle 
arrive  at  the  central  tube  as  weak  vortices.  Consequently, 
each  tube  behaves  more  or  less  as  if  it  were  independently 
subjected  to  a  turbulent  harmonic  flow.  As  K  increases,  not 
only  the  turbulence  level  but  also  the  interaction  between 
the  wakes  of  the  various  cylinders  increases.  There  is  a 
certain  amplitude  of  oscillation  beyond  which  neither  the 
interaction  of  the  wakes  nor  the  increase  of  the  turbulence 
level  affects  the  overall  force  acting  on  the  bundle.  A  com¬ 
parison  of  the  total  drag  force  acting  on  the  bundle  with  the 
sum  of  the  drag  forces  acting  on  each  cylinder  in  isolati  -n 
in  harmonic  flow  (at  the  corresponding  K  and  Re  values 
appropriate  to  each  tube)  shows  that  the  former  is  about  10 
percent  smaller. 

The  lift  coefficients,  as  given  by  Eqn.  (5),  for  the  two 
riser  configurations  are  shown  in  Figures  7c  and  7d.  Apparently, 
both  configurations  yield  nearly  identical  lift  coefficients. 

As  will  be  seen  later  the  rapid  decrease  of  the  lift  coefficient 
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with  increasing  K  is  in  conformity  with  that  for  a  single 
cylinder.  As  the  Keulegan-Carpenter  number  increases,  each 
member  of  the  bundle  sheds  an  increasing  number  of  vortices 
with  different  strengths  and  frequencies. 

Furthermore  the  shielding  or  proximity  effect  of  the  tubes 
results  in  a  highly  confused  vortex  system.  This  in  turn 
decreases  the  coherence  of  the  regular  vortex  shedding,  re¬ 
sulting  in  very  low  lift  coefficients.  Evidently,  a  greater 
understanding  of  the  lift  force  requires  the  determination  of 
the  lift  force  acting  on  each  member  individually  and  the 
consequence  of  tube  proximity. 


B.  DRAG  AND  INERTIA  COEFFICIENTS  FOR  THE  TWO  CYLINDER 
ARRANGEMENT 

The  drag  and  inertia  coefficients  for  N=1.5,  N=2,  N=2.5, 
and  N=3 . 5  (only  for  ci=0)  and  for  a=0°,  a=30°,  a=60°,  and  a=90° 
are  shown  in  Figures  8  through  20.  Each  figure  represents 
either  the  drag  or  the  inertia  coefficient  for  two  different 


flow  directions.  For  example.  Figure  8a  shows  the  drag 
coefficients  C^  and  C ^  for  ct=0°  and  N=1.5.  The  coeffi¬ 
cients  C^  is  for  the  flow  directed  from  passive  towards 
active  cylinder.  The  coefficient  C^  is  for  the  flow  dir¬ 
ected  from  the  active  towards  the  passive  cylinder.  In 
addition,  included  in  each  figure  is  the  drag  or  inertia 
coefficient  for  the  single  cylinder  of  the  same  size. 

Figure  8b  shows  the  inertia  coefficients  C^^  and  C  2  f°r 

a=0°  and  N=1.5.  The  coefficients  C  .  and  C  _  are  related  to 

ml  m2 

the  flow  direction  in  the  same  manner  as  the  coefficients 


Cdl  and  Cd2- 
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Figures  8a,  9a,  10a,  and  11a  show  that  for  the  tandem 
arrangement  corresponding  to  a=0°  the  drag  coefficient 
is  always  larger  than  .  The  difference  between  the  two 
drag  coefficients  decreases  with  increasing  N  (larger  spacing) . 
This  finding  is  in  conformity  with  the  measurements  reported 
by  Zdravkovich  /— 5 J  for  the  tandem  arrangement  in  steady 
flow.  For  relatively  small  spacings  (smaller  N)  the  wake  of 
the  passive  cylinder  is  highly  turbulent.  This  causes  earlier 
separation  on  the  active  cylinder  and  hence  in  smaller  C^ 
values.  One  may  also  express  the  same  fact  by  stating  that 
the  active  cylinder  is  shielded  by  the  passive  cylinder  when 
the  flow  sees  the  passive  cylinder  first.  As  the  spacing 
between  the  two  cylinders  increases,  and  gradually 

approach  each  other  and  single  cylinder  value.  Even  for 
N=3 . 5 ,  C^2  is  slightly  larger  than  because  of  the  afore¬ 

mentioned  effect  of  turbulence  in  the  wake  of  passive  cylinder. 

The  inertia  coefficients  C  ^  and  for  the  tandem  arrange¬ 

ment  are  shown  in  Figures  8b,  9b,  10b,  and  lib.  Unlike  the 
drag  coefficient,  the  inertia  coefficient  C  ^  is  always 
smaller  than  C  .  This  is  not  entirely  unexpected.  All  drag 
and  inertia  evaluations  for  single  cylinders  in  waves  and 
harmonic  flows  have  shown  that  /— 20_/  the  drag  and  inertia 
coefficients  follow  opposite  trends.  In  other  words,  when 
the  drag  coefficient  increases  the  inertia  coefficient  de¬ 
creases  and  vice  versa.  The  reason  for  this  is  embedded  in 


the  formulation  of  the  Morison's  equation  (See  Eq.  (1)) .  When 
the  phase  angle  between  the  velocity  and  the  maximum  force 
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is  small  the  tctal  force  is  primarily  due  to  the  velocity- 
square  dependent  drag  ana  the  contribution  of  the  inertial 
force  is  negligible.  In  fact  such  flows  are  termed  drag- 
dominated  flows.  When  the  phase  angle  between  the  velocity 
and  maximum  force  is  large,  the  total  force  is  primarily  due 
to  the  acceleration-dependent  inertia.  Such  flows  are  termed 
inertia-dominated  flows.  Thus,  it  is  apparent  that  the  drag 
and  inertia  coefficients  should  follow  opposite  trends. 

As  the  relative  spacing  increases  (See  Fig.  lib)  C  ^  and 
approach  each  other.  However,  the  difference  between  two 
coefficients  remains  relatively  larger  than  that  between 
and  C^2 •  The  experimental  fact  is  that  the  spacing  between 
the  two  cylinders  in  tandem  arrangements  must  be  considerably 
larger  than  3.5  diameters  for  the  inertia  coefficients  to 
become  independent  of  the  proximity  effects.  It  is  also  clear 
that  structures  in  close  proximity  can  cause  relatively  more 
significant  changes  in  the  inertia  coefficient  than  in  the 
drag  coefficient. 

Figures  8b,  9b,  10b,  and  lib  also  show  that  the  inertia 
coefficient  for  single  cylinder  (dashed  lines)  is  nearly  the 
same  as  C  2 •  ln  other  words,  when  the  flow  sees  the  active 
cylinder  first  the  proximity  effect  on  the  inertia  coefficient 
is  negligible.  It  is  also  apparent  that  when  the  flow 
sees  the  passive  cylinder  first  the  inertia  coefficient  C  ^ 
is  considerably  increased.  The  important  conclusion  to  be 
reached  on  the  basis  of  these  results  is  that  the  inertia 
coefficient  of  a  shielded  cylinder  may  increase  considerably 
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and  the  drag  coefficient  may  decrease  significantly. 

The  other  extreme  of  the  two  cylinder  arrangement  is  the 
side-by-side  arrangement.  The  discussion  of  this  case  will 
be  taken  up  next  since  the  force  coefficients  for  all  other 
intermediate  values  of  a  follow  trends  intermediate  to  the 
tandem  and  side-by-side  arrangements. 

Figures  18a,  19a,  and  20a  show  and  C^2  for  N=1.5,  N=2, 
N=2 . 5  for  the  side-by-side  arrangement.  Several  facts  are 
immediately  apparent  from  these  figures.  For  spacings  as  small 
as  1.5  diameter  the  drag  coefficients  are  not  greatly  in- 
cluenced  by  the  proximity  effects.  Secondly,  and  C^2 

are  very  close  to  that  for  the  single  cylinder.  The  slight 
difference  in  the  lower  range  of  Keulegan-Carpenter  numbers 
may  be  attributed  to  the  increased  blockage  in  the  tunnel  and 
experimental  uncertainty  which  was  found  to  be  about  5  per¬ 
cent.  Experiments  /  5_/  with  steady  flows  about  the  same 
type  of  two-cylinder  arrangement  have  also  shown  that  the 
proximity  effects  are  negligible  for  spacings  larger  than 
about  1.5  diameter. 


The  inertia  coefficients  and  Cm2  shown  in  Figures  18b, 

19b,  20b  exhibit  a  somewhat  different  behavior.  Even  though 
they  are  nearly  identical  for  any  given  spacing  they  are  con¬ 
siderably  larger  than  that  for  the  single  cylinder.  Only  for 
N=2.5  that  they  approach  the  single  cylinder  values. 

The  two  other  arrangements  intermediate  to  those  discussed 


above  behave  as  would  be  expected.  Both  the  drag  and  inertia 
coefficients  shown  in  Figures  12  through  17  vary  between  the 
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cases  for  ct=0°  and  ct=90°  and  both  gradually  approach  the 
corresponding  single  cylinder  values.  The  conclusions  re¬ 
garding  the  effect  of  shielding  remain  valid  for  the  inter¬ 
mediate  arrangements. 


C.  LIFT  COEFFICIENTS  FOR  THE  TWO  CYLINDER  ARRANGEMENT 

As  in  the  case  of  drag  and  inertia  coefficients,  first  the 
tandem  arrangement,  then  the  side-by-side  arrangement,  and 
finally,  the  intermediate  arrangements  will  be  discussed. 

For  the  tandem  arrangement  cjji=Cjj2  an<^  <“l3=CL4  ^ee 
6a  and  6b)  .  The  coefficients  C^  or  C^  are  shown  in  Figures 
21-24  and  the  coefficients  or  CL4  are  shown  in  Figures 

25-28.  Also  shown  in  each  figure  is  the  lift  coefficient  for 
a  single  cylinder  of  the  same  size. 

The  data  showed  that  the  lift  coefficient  is  strongly 
affected  by  the  cylinder  proximity,  irrespective  of  the 
direction  of  the  flow.  For  Keulegan-Carpenter  numbers  larger 
than  about  20  the  lift  coefficients  C^  and  C  ^  remain  prac¬ 
tically  constant.  Only  for  very  small  Keulegan-Carpenter 
numbers  that  C  ^  and  increase  rapidly,  in  conformity  with 

that  for  a  single  cylinder.  The  reason  for  this  is  that  for 
Keulegan-Carpenter  numbers  near  10 ,  at  most  two  vortices  are 
shed  during  a  given  half-cycle.  These  vortices  do  not  travel 
1.6  diameters  in  one  direction  as  seen  from  the  simple  relation 
that  K=2ta/D  where  A  is  the  amplitude  of  the  flow  oscillation. 
Consequently,  the  vortices  shed  from  the  active  cylinder  are 
only  partially  blocked  by  the  passive  cylinder  for  K  smaller 
than  about  10.  Nevertheless,  the  effect  of  cylinder  proximity 
is  ever  present  as  seen  from  a  comparison  of  the  data  with 


that  for  the  single  cylinder.  As  the  spacing  increases,  the 
lift  coefficients  gradually  approach  that  of  the  single  cylin¬ 
der  as  seen  in  the  Figures  24  and  28. 

In  the  case  of  side-by-side  arrangement  CL]_=CL3  and 

CL2=CL4*  coe^^^-c^ents  or  Cj.3  are  s^own  ^-n  Figures  29- 

31  and  the  coefficients  C L2  or  CL4  /  in  Figures  32-34.  The 
lift  coefficients  for  this  particular  arrangement  have  been 
measured  for  only  N=1.5,  N=2,  and  N=2.5. 

As  the  flow  oscillates  about  the  cylinders  arranged  in 
side-by-side  position,  the  blockage  introduced  by  the  proximity 
of  the  cylinders  forces  the  fluid  to  flow  above  the  top  and 
below  the  bottom  cylinder.  In  other  words,  the  velocity  on 
the  extreme  sides  of  the  passive  and  active  cylinder  increases 
while  the  velocity  between  the  two  cylinders  decreases.  This 
in  turn  increases  the  pressure  on  the  sides  of  the  cylinders 
nearest  each  other  and  decreases  the  pressure  on  the  opposite 
sides.  Consequently,  or  should  be  smaller  than  C^2 

or  C^.  A  comparison  of  the  Figures  29-31  with  Figures 
32-34  shows  that  C  or  is  in  fact,  always  smaller  than 

CL2  or  Cl4 •  As  the  relative  spacing  increases  all  four 
coefficients  gradually  approach  that  for  a  single  cylinder  as 
seen  from  Figures  31  and  34.  These  results  are  in  conformity 
with  those  obtained  with  steady  flow  about  a  side-by-side 
arrangement  of  circular  cylinders.  /-5_7 

For  arbitrary  arrangements  one  must  deal  with  four  dif¬ 
ferent  lift  coefficients,  namely,  ,  C^,  ,  and  CL4 . 

The  coefficient  is  shown  in  Figures  35-37  for  a=30°  and 
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for  N=1 . 5 ,  N=2 ,  and  N=2.5.  Similarly;  C 


is  shown  Figures 


L2 

38-40;  C.^ >  in  Figures  41-43;  and  finally,  ,  in  Figures 
44-46.  The  corresponding  lift  coefficients  for  a=60°  arrange¬ 
ment  is  shown  in  Figures  47-58. 

A  comparison  of  Figures  35-37  with  Figures  38-40  shows 
that  is  always  larger  than  •  The  reason  for  this  is 

that  the  vortices  shed  from  the  passive  cylinder  pass  over 
the  top  of  the  active  cylinder  and  thereby  reduce  the  pressure 
at  the  top  of  the  active  cylinder.  The  coefficients  and 

0^2  gradually  approach  that  for  a  single  cylinder  as  the 
spacing  increases. 

When  the  flow  first  sees  the  active  cylinder,  C  ,  remains 

L  J 

relatively  larger  than  because  of  the  presence  of  the 

vortices  near  the  top  of  the  active  cylinder  which  were  shed 
during  the  previous  half  cycle.  As  the  spacing  between  the 
two  cylinders  increases  both  and  rapidly  approach  the 

lift  coefficient  for  a  single  cylinder. 

For  ct=60°  arrangement  all  four  lift  coefficients  are 
practically  unaffected  by  the  cylinder  proximity  even  for  a 
spacing  as  small  as  1.5  diameters  and  they  are  nearly  equal 
to  that  of  an  isolated  cylinder.  Evidently,  the  effect  of 
cylinder-cylinder  proximity  on  the  lift  coefficient  for  a 
harmonically  oscillating  flow  is  confined  to  relatively  small 
spacings  and  angular  positions. 

The  frequency  of  the  lift  force  oscillations  has  been 
expressed  in  terms  of  the  Strouhal  number  and  evaluated  for 
all  cylinder  arrangements.  The  results  have  shown  that  the 
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Strouhal  number  remains  constant  at  an  average  value  cf  about 
0.2.  The  only  exception  to  that  was  the  case  corresponding 
to  a=30°  and  N=1.5. 

Figures  59  and  60  show  the  Strouhal  number  as  a  function 
of  K  for  a=0°  and  N=1 . 5 .  As  N  increased  the  Strouhal  numbers 
for  the  two  flow  directions  rapidly  approached  each  other  and 
a  mean  value  of  0.2. 

Figures  61  and  62  show  the  Strouhal  number  for  a=90°  and 
N=2 .  Once  again  the  lift  coefficients  for  the  two  directions 
of  flow  are  identical  and  gradually  approach  a  value  of  about 
0.2.  In  the  case  of  ct=30°,  the  wakes  interfere  in  an  asymme¬ 
tric  manner  and  the  Strouhal  number  shows  a  strong  dependence 
on  the  Keulegan-Carpenter  number.  Figures  63-68  show  the 
Strouhal  number  as  a  function  of  K  for  various  values  of  N 
for  a=30°.  Apparently,  the  Strouhal  number  decreases  from 
values  as  large  as  0.4  to  about  0.15  as  K  increases.  For 
spacings  larger  than  about  two  diameters  Strouhal  number 
decreases  to  about  0.2. 
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VI .  CONCLUS IONS 


An  extensive  investigation  of  the  forces  acting  on  tube 
bundles  and  arrangements  of  two  cylinders  has  wa  anted  the 
following  conclusions: 

1.  Experiments  with  two  particular  riser  configurations 
have  shown  that  the  inertia  coefficient  is  considerably 
larger  than  that  predicted  by  the  potential  theory 

and  shows  that  some  fluid  mass  is  entrapped  within  the 
bundle  as  a  consequence  of  the  solidification  of  the 
tube  configuration. 

2.  For  tandem  arrangement  of  two  cylinders  the  drag  and 
inertia  coefficients  depend  on  both  the  Keulegan- 
Carpenter  number  and  the  relative  spacing  of  the 
cylinders.  As  the  amplitude  of  flow  oscillation  be¬ 
comes  comparable  or  smaller  than  the  gap  between  the 
two  cylinuers,  the  drag  and  inertia  coefficients 
gradually  approach  those  corresponding  to  an  isolated 
cylinder,  i.e.,  the  wakes  of  the  two  cylinders  do  not 
interfere  with  each  other. 

3.  The  drag  and  inertia  coefficients  for  the  side-by-side 
arrangement  exhibit  a  similar  behavior.  For  a  relative 
spacing  larger  than  about  2.5  diameters,  the  cylinders 
behave  as  if  they  were  independent. 

4.  The  quantification  of  the  lift  force  for  an  arbitrary 
arrangement  requires  four  different  lift  coefficients, 
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depending  on  the  direction  of  the  flow  and  the  direc¬ 
tion  of  the  force.  These  coefficients  have  been  eval¬ 
uated  and  it  has  been  shown  that  the  effect  of  cylinder- 
cylinder  proximity  is  confined  to  relatively  small 
spacings  and  angular  positions. 

5.  The  Strouhal  number  for  practically  all  two-cylinder 
arrangements  remained  nearly  constant  at  about  0.2. 

Only  in  the  case  of  30°  orientation  did  the  strong 
asymmetric  interference  of  the  cylinder  wakes  result 
in  a  Keulegan-Carpenter-number  dependent  Strouhal 
number,  ranging  from  0.4  to  0.2.  However,  as  the 
spacing  increased  the  Strouhal  number  again  approached 
a  mean  value  of  0.2. 
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Fig.  7a  C  and  C.  versus  K  for  riser 
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Fig.  7c  C  vs  K  for  riser 
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Fig.  8a  C.-.  and  versus  K  for  N=1.5  and  a=0° 


Fig.  10a  C  .  and  C,9  versus  K  for  N=2.5  and 


Fig.  14a  C  ,  and  C,2  versus  K  for  N=2.5  and  a=30 


Fig.  17a  C.  and  C,_  versus  K  for  N=2.5  and  a=60 
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Fig.  18a  C,-.  and  CH~  versus  K  for  N=1.5  and  a=90 


Fig.  18b  C,.  and  versus  K  for  N=1.5  and  oc=90 


67 


Pig.  19b  C  -i  and  C  »  versus  K  for  N=2  and  ot=90 

J  ml  m  ) 


I 


Fig.  22  C  or  C  ,  vs  K  for  N=2  and 
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Fig.  23  C  .  or  C  0  vs  K  for  N=2.5  and 


Fig.  24  C  ,  or  C  2  vs  K  for  N=3.5  and 


Single  Cylinder 
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Fig.  25  C  ,  or  C  .  vs  K  for  N=1.5  and 
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Fig.  26  C  _  or  C  .  vs  K  for  N=2  and 


Fig.  28  C  ^  or  C  .  vs  K  for  N=3.5  and 


Single  Cylinder 
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Fig.  29  CT ,  or  CT,  vs  K  for  N=1.5  and  a=90 
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Fig.  30  CT or  CT vs  K  for  N=2  and  a=90 
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or  CT  ->  vs  K  for  N=2.5  and  a-90 
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Fig.  32  C  9  or  C  .  vs  K  for  N=1.5  and  a=90 
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Fig.  33  CT 0  or  C  .  vs  K  for  N=2  and  a=90 


Fig.  34  C  ~  or  C  .  vs  K  for  N=2.5  and  a=90 
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Fig.  36  C  ,  vs  K  for  N=2  a=30 


Fig.  37  C  .  vs  K  for  N=2.5  and  a=30 
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Fig.  37  C  0  vs  K  for  N=1.5  and 
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Fig.  39  C  _  vs  K  for  N=2  and 
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Fig.  41  C.  vs  K  for  N=1.5  and  a=30 


92 


Fig.  42  C  _  vs  K  for  N=2  and  a=30 
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vs  K  for  N=2.5  and 


94 


Fig.  44  C  .  vs  K  for  N=1.5  and  cx=30 


vs  K  for  N=2  and 
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Fig.  46  C  .  vs  K  for  N=2.5  and 


Fig.  48  C  .  vs  K  for  N=2  and  a=60 


Fig.  49  C  vs  K  for  N=2.5  and  a=60 


Fig.  51  C  „  vs  K  for  N=2  and  a=60 
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Fig.  52  C  ~  vs  K  for  N=2  and  a=60 


Passive  Cyl 


Fig.  53  CT  i  vs  K  for  N=1.5  and  a=60 


Fig.  55  C  vs  K  for  N=2.5  and  a=60 


Fig.  54  C  _  vs  K  for  N=2  and  a=60 


Fig.  56  C  .  vs  K  for  N=1.5  and  a=60 
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Fig.  57  C  .  vs  K  for  N=2  and 


Single  Cylinder 


Fig.  58  C  vs  K  for  N=2.5  and  ct=60 


Fig.  60  St ^  and  St  vs  K  for  N=1.5  and 
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for  N=1 . 5  and  a=30 
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Fig.  67  St.  and  St,,  vs  K  for  N=2.5  and  a=30 
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vs  K  for  N=2 . 5  and  a=30 
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